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Identifying low-dimensional dynamics in Type-I edge-localised-mode 
processes in JET plasmas 
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Edge localised mode (ELM) measurements from reproducibly similar plasmas in the Joint European Torus 
(JET) tokamak, which differ only in their gas puffing rate, are analysed in terms of the pattern in the sequence 
of intcr-ELM time intervals. It is found that the category of ELM defined empirically as Type I - typically 
more regular, less frequent, and having larger amplitude than other ELM types - embraces substantially 
different ELMing processes. By quantifying the structure in the sequence of intcr-ELM time intervals using 
delay time plots, we reveal transitions between distinct phase space dynamics, implying transitions between 
distinct underlying physical processes. The control parameter for these transitions between these different 
ELMing processes is the gas puffing rate. 

Keywords: tokamak edge plasmas, edge localised modes, nonlinear time series, delay time plots 
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I. INTRODUCTION 

Enhanced confinement operating regimes (H-mode) in 
magnetically confined plasmas are accompanied by pulses 
of energy and particle release known as edge localised 
modes (ELMs)i~— . At steady state, a magnetically con- 
fined tokamak plasma comprises a family of nested mag- 
netic flux surfaces in a smooth, or laminar state. ELM- 
ing constitutes a relaxation process, for the edge region 
of H-mode plasmas, which encompasses an initial trig- 
ger for linear MHD instability evolving into a fully non- 
linear detached state, such that structures propagate to 
the first wall where they generate recombination radia- 
tion. In parallel, local temperature and pressure gra- 
dients evolve rapidly. The onset of ELMing accompa- 
nies a sharp transition in the global state of the toka- 
mak plasma, and changes in observed ELM character 
reflect changes in externally applied drive such as gas 
puffing and heating. Control, mitigation and prediction 
of the occurrence of large Type I ELMs are central cha- 
llenges for magnetic confinement fusion plasma physics. 
There are many active experimental campaigns in this 
area- - — , particularly in support of the future ITER toka- 
mak, for which the consequences of uncontrolled Type 
I ELMs may be unacceptable 4 ^. While successful theo- 
ries for some component elements of the ELMing pro- 
cess have been constructed, there is currently no com- 
prehensive first principles model that incorporates all 
of the physical effects that are known to contribute to 
the ELMing process. ELM categorisation is primarily 
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phenomcnological 3 - - — , furthermore it is not always easy 
to discriminate in real time between Type I and, say, 
Type III ELMing. Hitherto only a few paper s 10 ' 11 have 
addressed measured ELM sequences as the pulsed out- 
puts of a nonlinear system, a field where generic analy- 
sis techniques are well developed and potential links to 
ELMing have long been apparent—. Characterisation of 
ELMing processes by applying dynamical systems theo- 
ry to the data offers a fresh avenue to understanding, 
prediction and control, and may help identify some of 
the key properties that models for Type I ELMing must 
embody. Here we take the first steps. 

D. Ruelle and F. TakensH initiated a classical sce- 
nario for the transition from ordered to disordered flow 
in fluids with increasing driving control paramete r 14 ' 15 . 
This has been observed in Rayleigh-Benard convection 
in fluids^—, and in drift wave turbulence 2 ^ and flute 
instabilities in plasmas^ 2 -. Oscillatory behaviour arises 
cither if there is a constant of the motion, or if there is 
a limit cycle onto which the system dynamics is attrac- 
ted in the presence of damping or dissipation. In the 
present case, where the system is the plasma undergoing 
the ELMing process, the nature and number of the re- 
levant phase space co-ordinates is not known from first 
principles. Progress towards their identification can ne- 
vertheless be made by applying techniques of dynamical 
systems analysis to visualize changes in the topology of 
the phase space. A convenient method is that of 'de- 
lay plots', that is, to plot the successive time intervals 
between crossings of a surface of section in the phase 
spaced— . 

In this article we report the application of delay plots 
to the measured time intervals or waiting times between 
successive ELMs. We consider ELM sequences from 
six similar plasmas in the JET tokamak, including JET 
plasma 57865 where the H-mode closely approaches an 



ITER operating regime with respect to some, but not 
all, key dimensionless parameters^. We obtain evidence 
that Type I ELMing in these plasmas exhibits transi- 
tions between processes with distinct physical analogues, 
dependent on the value of the gas puffing rate as con- 
trol parameter. In all six plasmas the toroidal magnetic 
field density is 2.7T, the plasma current is 2.5MA, neu- 
tral beam and ion cyclotron resonance heating power are 
13.5MW and 2.0MW respectively, and the H g $ confine- 
ment factor is in the range 0.87 to 1.0. In all six plasmas, 
gas puffing terminates at 23.3s and neutral beam heat- 
ing is ramped down from 23.5 to 24.5s. The differences in 
Type I ELM character are largely determined by the dif- 
ferent levels of externally applied gas puffing. The inten- 
sity of the D a signal, which sometimes saturates, is not 
necessarily a reliable proxy for the magnitude of the un- 
derlying ELM plasma phenomenon, whereas occurrence 
times are well defined. ELM occurrence and ELM wait- 
ing times are the primary physical indicators addressed 
in the present study. The moment of occurrence of each 
ELM is inferred from the D a datasets using an algorithm 
similar to that described ir*ii, which exploits the steep 
leading edge of each ELM. This procedure generates a 
sequence of event times t n for each nth ELM, and hence 
inter-event times St n = t n — t n -i- These sequences are 
used to construct delay plots, which are known22r— to 
capture aspects of the topology of the unknown underly- 
ing phase space evolution of the system. 



II. RESULTS AND DISCUSSIONS 

Figures 1 and 2 show measured Type I ELM signals for 
a sequence of six JET H-modc plasmas 578nrn, where 
nm is 72, 71, 70, 65, 67, and 69 in order of increasing 
magnitude and duration of the gas puffing rate, shown in 
Fig. 3, which is the key external control parameter. The 
upper trace in each panel of Figs. 1 and 2 plots the 
time-evolving intensity of Lyman alpha recombination 
radiation from deuterium, D a , measured by a camera 
directed at the inner divertor, normalised by the mean 
measured intensity. The two groupings of three plasmas 
are at lower (Fig. 1) and higher (Fig. 2) gas puffing 
rates. At lower gas puffing rates (Fig. 1) the ELM sig- 
nal intensity is roughly the same across each time series, 
whereas at higher gas puffing rates (Fig. 2) this shows 
a rich structure. We will investigate this structure by 
sorting the ELM events that are used to construct the 
time series of inter-ELM time intervals (delay times), in 
terms of whether they exceed a threshold in signal in- 
tensity; the thresholds used arc indicated by horizontal 
lines on the ELM time series (top panel in Figs 1 and 
2). Each nth Type I ELM that has signal intensity ex- 
ceeding a given threshold then forms a set of events at 
time t„ with the delay between events St n = t n — t n —\. 
The middle panels of Figs. 1 and 2 show the delay plots 
for a given threshold, that is, 5t n+ i versus St n . The D a 
signal intensity for the ELM at t n is indicated by colour 



coding. These delay plots reflect the topology of the sys- 
tem phase space. For a trajectory that is approximately 
singly periodic the delay plot will exhibit a single con- 
centration of points on the 5t n +i = 5t n line, centred 
on the mean period r = 5t n +i = St n . The spread of 
points about the mean period reflects a combination, in 
unknown proportions, of intrinsic and extrinsic sources 
of irregularity in a quasi-regular process, and determines 
the practical resolution limit of this method. A period- 
two oscillation will generate two concentrations of points, 
symmetrically placed either side of the 5t n +\ = St n line. 
Dynamical switching between one period n and another 
at Ti will generate four concentrations of points: at the 
two distinct periods T\ and ti on the bt n -\-i = St n line, 
and at two locations symmetrically placed either side of 
the line, at (8t n -\-i,5t n ) coordinates (ti,T2) and (t2,Ti). 
The number of ELMs evaluated in these six JET plas- 
mas ranges between 79 and 197. The mean inter-ELM 
time interval is in the range 25 to 60ms. The delay plots 
in Fig. 1 are insensitive to the threshold, in marked con- 
trast to Fig. 2, suggesting that these reflect distinct pro- 
cesses. In Fig. 1, plasmas with successively greater gas 
puffing rates are shown from left to right. We can see 
that increased gas puffing causes the ELMing process to 
bifurcate from singly periodic (57872), via transitional 
behaviour (57871), to a situation where two periods are 
present (57870) together, with the plasma switching be- 
tween them. This behaviour is approximately analogous 
to that of small amplitude oscillations of two weakly cou- 
pled pendulums with different natural frequencies. It is 
also apparent that a longer delay time 8t„ before an ELM 
correlates statistically with a larger D a signal intensity. 
The bottom pair of plots in each panel of Figs. 1 and 2 
displays the probability density functions (pdfs) for the 
distributions of measured St n for the ELM time series vi- 
sing the same amplitude thresholds as for the delay plots; 
in Fig. 1, unlike Fig. 2, these two panels are identical. 

We now turn to Fig. 2 which corresponds to higher 
overall levels of gas puffing rate. It displays a transition 
in the ELMing process as the gas puffing rate is increased, 
which is different to that seen in Fig. 1. Each ELM 
with large D a signal intensity is statistically likely to 
be rapidly followed by a population of postcursor ELMs 
with smaller D a signal intensity. The likelihood of a 
postcursor ELM, and their number, increases with gas 
puffing rate. As a consequence, the delay plots construc- 
ted for different thresholds now, unlike Fig. 1, show dif- 
ferent structure. At relatively low gas puffing rate (left 
hand plots) most delays fall within a single group on 
the St n+ i = St n line. However when the threshold is 
reduced, smaller postcursor events begin to feature in 
the time series of delays and result in populations (lines 
parallel to the axes) far from the 5t n +i = St n line, and 
a new, narrowly constrained group on the 5t n +i = St n 
line at small (<5i n +i, 5t n ). As the gas puffing rate is in- 
creased, these small postcursor events come to dominate 
numerically. It is noteworthy that whereas ELMs with 
large signal amplitude exhibit a broad inter-ELM time 
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Figure 1. ELM characteristics of three similar JET plasmas 57872, 57871, 57870 at lower gas puffing rates, showing for each 
plasma: (top of each panel) the time trace of D a signal intensity, displaying also the two amplitude thresholds used for the centre 
and bottom plots; (centre of each panel) delay plots for ELMs, with amplitude colour coded above the higher (lower) threshold 
on the left (right); (bottom of each panel) corresponding probability density functions for the distributions of measured 5t n for 
the ELM time series, using the same amplitude thresholds as for the delay plots; the red and blue curves represent different 
binning of the same data. The three plasmas are ordered, from the left, in terms of increasing magnitude of gas puffing, see 
Fig. 3. 




19 20 21 22 

t(s) 

0.15 




0.05 0.1 0.15 
St (s) 



0.05 0.1 0.15 
St (s) 





n 


Mil 


II] 


llll 


F 


'lasma 


57 

T 


867' 

fllll 


& 




m- 




li 


Mil 


111 


# 


I 




0.1 0.2 

St„(s) 



0.1 0.2 

St (s) 




0.2 0.4 

St (s) 



10 




0.2 0.4 

St (s) 



600 
f 400 

Ll_ 

£200 
Q 



0.1 0.2 

St <s) 









A 


0.005 0.01 
8t n (S) 



0.1 0.2 

St (s) 



Figure 2. As Fig. 1, for three similar JET plasmas 57865, 57867, 57869 at higher gas puffing rates. The three plasmas are 
ordered, from the left, in terms of increasing magnitude of gas puffing, see Fig. 3. The bottom panels from JET plasmas 57867 
and 57869 also include an inset panel displaying the sharp peak in the PDF. The population in this sharp peak increases with 
the gas puffing rate, and the average period r — 6.7 ± 6.6 x 10 -2 (ms). 



interval distribution, the distribution of the postcursors 
is very sharply defined and is invariant between the three 
JET plasmas, see Fig. 2 bottom panels. Its inverse de- 
fines a potentially important characteristic frequency of 
the ELMing process. This process, as seen in the delay 
plots, is analogous to random large amplitude transient 
impulses driving a system that has a narrowband reso- 



nant frequency response. 

Figure 3 displays the gas puffing rates for all six JET 
plasmas. The clear changes in ELMing displayed in Fig. 
1, and for JET plasmas 57867 to 57869, arise under com- 
paratively small changes in gas puffing, while there is a 
relative large step (a factor of approximately two) be- 
tween 57865 and 57867. Other ELM interval dynamics 
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be assisted by characterizing the measured sequences of 
inter-ELM time intervals in this way. 
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Figure 3. Time trace of gas puffing rate, F, in particles per 
second, which is the primary external control parameter for 
the six otherwise similar JET plasmas: ordered, from the 
bottom, in terms of increasing magnitude. 



are in principle possible for other gas fuelling rates, espe- 
cially for fuelling rates between those of 57865 and 57867, 
for these otherwise identical plasma operating regimes. 

Some previous experiments have observed that Type I 
ELM frequency (mean inter-ELM interval) increases with 
gas puffing rate^ 7 -. Moreover, early theoretical studies^ 
suggested that it might be possible to explain the ex- 
perimentally observed transition from Type I to Type 
III ELMy H-mode triggered by strong gas puffing, as 
well as the subsequent increase in ELM frequency and 
deterioration of plasma confinement, as a transition from 
second to first stability (either ideal or resistive modes). 
However, there is still no widely accepted model for the 
overall ELMing process or processes. 



III. SUMMARY 

We have exploited the similarity of these six JET plas- 
mas which all have exceptionally long duration ~ 5s of 
the quasi-stationary ELMing process, and which appear 
to have only one effective control parameter, the gas puff- 
ing rate. These particular experiments yield a sufficient 
number of ELMs and inter-ELM times, to enable us to 
apply the delay plot technique to characterize the dy- 
namics. 

There exists an increasing number of ELMing regimes. 
These are typically characterized phenomcnologically 
and in terms of bulk plasma parameters. We believe 
that ELM interval analysis of the kind presented here, if 
applied more widely, will help shed light on such transi- 
tions in confinement phenomenology in tokamak plasmas. 
In particular, demonstrating and quantifying the effec- 
tiveness of ELM control and mitigation techniques will 
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